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Abstract: 


Electrons  moving  in  graphene  behave  as  massless  Dirac  fermions,  and  they  exhibit 
fascinating  low-frequency  electrical  transport  phenomena.  Their  dynamic  response, 
however,  is  little  known  at  frequencies  above  one  terahertz  (THz).  Such  knowledge  is 
important  not  only  for  a  deeper  understanding  of  the  Dirac  electron  quantum  transport, 
but  also  for  graphene  applications  in  ultrahigh  speed  THz  electronics  and  IR 
optoelectronics.  In  this  paper,  we  report  the  first  measurement  of  high-frequency 
conductivity  of  graphene  from  THz  to  mid-IR  at  different  carrier  concentrations.  The 
conductivity  exhibits  Drude-like  frequency  dependence  and  increases  dramatically  at 
THz  frequencies,  but  its  absolute  strength  is  substantially  lower  than  theoretical 
predictions.  This  anomalous  reduction  of  free  electron  oscillator  strength  is  corroborated 
by  corresponding  changes  in  graphene  interband  transitions,  as  required  by  the  sum  rule. 
Our  surprising  observation  indicates  that  many-body  effects  and  Dirac  fermion-impurity 
interactions  beyond  current  transport  theories  are  important  for  Dirac  fermion  electrical 
response  in  graphene. 
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Graphene  provides  a  unique  material  system  to  study  Dirac  fermion  physics  in  two 
dimensions.  Researchers  have  demonstrated  in  graphene  exotic  Dirac  fermion 
phenomena  ranging  from  anomalous  quantum  Hall  effects  ’  to  Klein  tunneling  in  low 
frequency  (DC)  electrical  transport.  They  also  observed  an  optical  conductance  defined 
by  the  fine  structure  constant4’5  and  gate-tunable  infrared  absorption  in  Dirac  fermion 
interband  transitions  ’  .  Situated  between  low-frequency  electrical  transport  and  interband 
optical  excitation  is  the  spectral  range  dominated  by  free  carrier  intraband  transitions, 
studies  of  which  can  provide  new  insight  into  both  Dirac  fermion  transport  and  interband 
transitions.  This  free  carrier  dynamics  response  is  also  expected  to  play  a  key  role  in 
future  development  of  ultrahigh-speed  electronics  at  THz  frequencies  and  THz-to-mid-IR 
optoelectronic  devices8,9.  However,  despite  its  importance,  experimental  study  of  free 
carrier  response  of  graphene  in  the  THz  to  mid-IR  spectral  range  has  been  lacking. 

Quantum  theories  of  electron  dynamic  response  in  graphene  have  been  developed  by 
several  groups  10'15.  They  predict  that,  within  the  framework  of  Boltzmann  transport 
theories,  high  frequency  conductivity  of  monolayer  graphene  has  a  Drude 

form,  a (a>)  = - — - ,  where  co  is  the  frequency  and  T  the  scattering  rate.  The 

7r((o  +  iT) 

y  ^  I — j — - 

prefactor  D,  known  as  the  Drude  weight,  has  the  value  D  =  — — J/nn  when  electron- 

h 

electron  interactions  are  neglected,  with  vF  =  l.lxl06m/,v  being  the  Fermi  velocity2  and 
Ini  the  carrier  density.  This  is  distinctly  different  from  classical  materials  where 
7rn  • 

D  = -  .  The  predictions  on  graphene  Dirac  fermion  dynamics  have  never  been 

m* 

directly  tested  experimentally.  In  fact,  it  was  suggested  from  measurements  of  interband 
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f. : 

transitions  that  D  could  deviate  from  the  predicted  value.  There  are  also  theoretical 
studies  showing  that  electron-impurity  interactions  could  result  in  carrier  dynamics  in 
graphene  beyond  the  Boltzmann  description  16’17,  and  that  inclusion  of  electron-electron 
interactions  could  dramatically  reduce  the  Drude  weight  D  18. 

Here,  we  report  spectroscopic  measurements  of  graphene  conductivity  from  terahertz 
to  mid-IR  for  different  electron  and  hole  concentrations.  We  confirm  that  the 
conductivity  from  free-carrier  response  indeed  has  a  Drude-like  frequency  dependence. 
In  addition,  we  have  been  able  to  determine  the  Drude  weight  D  and  scattering  rate  T  of 
doped  graphene  directly  and  separately  for  the  first  time.  We  observe  an  electron-hole 
asymmetry  for  D,  and  find  that  values  of  D  is  substantially  lower  than  the  theoretical 
^2 

prediction  D  =  — — Jx \n\  10"12.  It  indicates  that  a  better  understanding  of  electric 
h 

transport  and  Drude  conductivity  need  to  include  electron-electron  and  electron-impurity 
interactions  beyond  the  current  theories  for  carrier  dynamics.  This  reduction  of  Drude 
weight  (arising  from  free  carrier  intraband  transitions)  is  connected  to  corresponding 
changes  in  interband  transitions  (6),  and  we  demonstrate  that  the  sum  rule  requiring  the 
integrated  intraband  (Drude)  conductivity  and  interband  conductivity  to  be  a  constant  is 
well  obeyed. 

In  our  study  we  used  Fourier  transform  infrared  spectroscopy  (FTIR)  to  measure  the 
transmission  spectra  of  graphene  samples  over  the  range  from  30  to  6000  cm'1  and 
deduce  the  frequency-dependent  conductivity  from  the  spectra.  Previously,  infrared 
spectroscopy  to  probe  graphene  monolayer  was  limited  to  >1000  cm'1  (wavelength  <10 
fj.m)  because  of  the  limited  size  of  exfoliated  graphene  4’6’7’19.  For  spectroscopy  at  longer 
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wavelength  up  to  -300  /rm ,  we  need  large  sample  size.  Here  we  used  large-area 


graphene  grown  by  chemical  vapor  deposition  (CVD).  (Fig.  la)  (see  methods) 

Fig.  lb  shows  the  DC  conductivity  as  a  function  of  the  gate  voltage  in  a  typical 
graphene  sample.  The  minimum  conductance  occurs  at  14  V,  which  defines  the  charge 
neutral  point  (CNP)  of  the  sample.  Carrier  density  in  graphene  is  related  to  the  gate 
voltage  by  n  =  7.5  x  1 010  (Vcnp  —  Vg )  cm~2V~]  2,  where  negative  (positive)  density 

corresponds  to  electron  (hole)  doping.  The  sample  had  an  initial  hole  doping  of  1.05x1 012 
cm'2  at  Vg  =  0,  and  a  hole  mobility  around  2700  cm2/V.s.  As  seen  from  Fig.  lb,  the  DC 
conductivity  for  hole  doping  is  reasonably  linear  with  Ini,  and  there  is  a  large  asymmetry 
between  electron  and  hole  conductivity  at  the  same  Ini.  In  previous  studies,  conclusions 
have  frequently  been  made  on  carrier  scattering  from  such  DC  conductivity  data  based  on 

c 7dc  =  D/ 7iY  and  the  assumption  that  D  =  — — Jn\n\  holds  exactly  as  theory  predicts21" 

h 

23.  However,  the  validity  of  this  approach  has  never  been  tested. 

An  independent  determination  of  D  and  Y  can  be  achieved  through  AC  conductivity 
measurements  using  IR  spectroscopy.  Fig.  lc  shows  a  difference  IR  absorption  spectrum 
of  hole-doped  graphene  (Vg  =  -70V)  in  reference  to  absorption  at  the  CNP  (Vg  =  14V). 
This  difference  absorption  spectrum  has  two  characteristic  features:  a  large  absorption 
increase  at  lower  wavenumbers,  which  reaches  over  15%  at  terahertz  frequencies  for  a 
graphene  monolayer,  and  an  absorption  reduction  over  a  broad  range  of  higher 
wavenumbers.  These  features  can  be  understood  qualitatively  from  the  gate-induced 
changes  in  intraband  and  interband  electronic  transitions  in  graphene  due  to  carrier 
doping  (Fig.  Id).  Free  carrier  conductivity  (from  intraband  transitions)  increases 
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dramatically  with  the  hole  doping.  It  peaks  at  zero  frequency  and  gives  rise  to  a 
substantial  absorption  increase  at  low  wavenumbers.  At  the  same  time,  however,  the 
interband  transitions  up  to  energy  2EF  (arrows  in  Fig.  Id)  become  forbidden  due  to  empty 
initial  states,  leading  to  a  reduction  of  absorption  in  the  broad  spectral  range  below  2EF 
(dashed  line  in  Fig.  lc). 

The  gate-induced  change  of  AC  conductivity  (referred  to  CNP  value).  Act  =  a- crcnp , 

can  be  obtained  readily  from  the  difference  IR  absorption  spectra  (see  supplementary 
information).  Fig.  2a  and  Fig.  2b  show  the  real  part  of  the  AC  conductivity  change  Act' 
in  the  low-wavenumber  range  (<450  cm'1)  for  different  electron  and  hole  concentrations. 
In  this  spectral  range  free  carrier  response  (i.e.,  an  increase  in  conductivity  with  carrier 
density)  dominates,  providing  direct  information  on  Dirac  fermion  electrical  transport. 
The  complete  set  of  Act'  spectra  in  Figs.  2a  and  2b  can  be  fit  (dashed  line)  by 
Acr(ct>)  =  iD  /  nico+iY)  —  iDcnp  /  Tiioj+iY  cnp )  using  the  Drude  form  for  o  and  ocnp  .  A 

finite  charge  neutral  point  response  ( Dcnp  and  T  )  accounts  for  inhomogeneous  electron 
and  hole  puddles  present  in  graphene.  In  the  fittings,  we  set  Ycnp  equal  to  that  when 
graphene  is  weakly  doped.  The  value  of  Dcnp  is  set  so  that  DCNP  /  7iTCNP  gives  the  DC 

conductance  at  the  charge  neutral  point.  From  the  Drude  fits,  we  obtain  directly  the 
scattering  rate  Y  and  Drude  weight  D  as  a  function  of  electron  and  hole  concentrations, 
shown  as  symbols  in  Fig.  2c  and  Fig.  2d,  respectively. 

From  the  values  of  D  and  T  extracted  from  the  IR  spectroscopy,  we  can  obtain 
ctdc(=D/ 7iY)  of  graphene  at  different  carrier  concentrations  (red  dots  in  Fig.  2e).  It 
agrees  nicely  with  the  directly  measured  DC  conductivity  in  our  sample  for  all  gate 
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voltages  (black  dots).  This  agreement  provides  a  consistency  check  of  the  Drude 
description. 

We  now  compare  our  experimentally  determined  T  and  D  to  theory  predictions. 
Figure  2c  shows  that  T  is  roughly  a  constant  at  different  hole  concentrations,  while  it  is 
higher  in  value  and  increases  with  the  carrier  concentration  for  electron  doping.  This 
electron/hole  concentration  dependence  of  T  is  not  predicted  by  theory,  and  cannot  be 
described  by  pure  unitary  scattering  or  charge  impurity  scattering  of  Dirac  fermions10'12’24. 
A  combination  of  different  scattering  mechanisms  seems  to  be  necessary  to  explain  the 
observation. 

While  theoretical  prediction  on  T  may  vary  depending  on  the  scattering  mechanisms 

^2 

assumed,  the  same  theoretical  value  for  the  Drude  weight  D  =  1W1153  is 

h 

always  obtained.  Comparison  of  theory  (blue  line)  and  experiment  (black  symbols)  in  Fig. 
2d,  however,  shows  an  appreciable  difference.  In  the  strongly  doped  region  (IVg-Vcnpl>20 
V)  where  the  Boltzmann  theory  is  supposed  to  be  most  reliable,  the  measured  D  is  lower 
than  the  theoretical  value  by  20-45%.  (The  experimental  uncertainty  is  less  than  10%.). 
An  electron  and  hole  asymmetry  is  also  present  in  the  measured  D.  The  experimentally 
observed  reduction  of  Drude  weight  is  quite  unexpected  and  it  suggests  that  the 
prevailing  graphene  transport  theory  is  not  complete. 

The  general  sum  rule  for  oscillator  strengths  in  solids  state  that  the  integrated 

/•oo 

conductivity  over  all  frequencies  is  a  constant,  i.e.  A (j\(o)d(o  =  0 .  In  the  case  of  non¬ 


interacting  free  carriers  described  by  prevailing  graphene  theories,  the  integrated 


intraband  absorption  increases  by  J  Acr'iK(ra  dco-  —  -  ^7r\n\  upon  carrier  doping 

2  2/j 


10-12 
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This  increase  is  compensated  exactly  by  an  integrated  interband  absorption 

\ A(J' inter  d(O=<J0'(2EF/h)  with  EF  =  kVF\[^H  and  ^0  =  ^  1 (Fig-  ld‘)  .  In 

our  case,  the  sum  rule  is  still  valid,  and  the  observed  anomalous  reduction  of  Drude 
weight  in  intra-band  transitions  should  be  accompanied  by  a  corresponding  change  in 
interband  transitions  at  large  wavenumbers.  Fig.  3a  and  3b  show  the  optical  conductivity 
difference  (referred  to  CNP)  at  different  hole  and  electron  concentrations  in  the  spectral 
range  600-6000  cm'1.  A  gated-induced  decrease  in  optical  conductivity  is  clearly  present 

2  E 

from  zero  to  a  cut-off  frequency  of — —  that  increases  with  carrier  doping,  as  the  theory 

h 

predicts.  However,  the  observed  interband  optical  conductivity  decrease  in  doped 
graphene  is  appreciably  less  than  cr0  (dashed  lines  in  Fig.  3  a  and  3b),  and  the  deviation  is 

larger  for  electron  doping  compared  with  hole  doping.  This  reduction  in  gate-induced 
absorption  decrease,  as  well  as  electron-hole  asymmetry  in  interband  transitions,  matches 
well  with  the  reduction  in  gate-induced  absorption  increase  and  electron-hole  asymmetry 
in  intraband  Drude  weight.  To  be  more  quantitative,  we  plot  in  Fig.  3c  the  absolute  values 
of  integrated  conductivity  change  for  interband  transitions  (black  symbol)  and  intraband 
transitions  (red  symbol,  which  is  proportional  to  the  Drude  weight  change)  at  different 
carrier  doping.  The  two  agree  almost  perfectly,  just  as  the  sum  rule  requires.  A  reduction 
of  gate-induced  interband  absorption  has  also  been  observed  previously  for  mechanically 
exfoliated  graphene6.  Based  on  the  sum  rule  that  we  established  experimentally,  it 
indicates  that  the  anomalous  Drude  weight  reduction  is  general  for  both  exfoliated  and 


CVD  samples. 


That  the  Drude  weight  D  is  significantly  smaller  than  the  predicted  value  of 
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D  =  — — \l7rn  ,  seen  directly  from  the  measurements  of  the  Drude  conductivity  and 
h 

confirmed  by  the  interband  absorption  spectra,  shows  that  modifications  on  the  prevailing 
transport  theory10'12  is  needed  to  describe  Dirac  fermions  in  graphene.  Inclusion  of 
electron-electron  interactions  and  improved  treatment  of  electron-impurity  interactions 
may  be  necessary.  It  has  already  been  shown  in  Ref.  18  that,  unlike  conventional  massive 
fermions,  the  Drude  weight  of  massless  Dirac  fermions  can  be  strongly  reduced  by 
electron-electron  interactions  due  to  pseudo-spin  physics,  although  this  theory  predicts  a 
much  larger  reduction  of  the  Drude  weight  (over  80%)  than  what  we  have  observed.  In 
addition,  the  electron-impurity  interactions  may  play  a  role  and  contribute  to  the  observed 
electron-hole  asymmetry17.  It  is  also  possible  that  reduction  of  D  is  a  consequence  of 
electron  or  hole  localization,  which  decreases  the  effective  ‘free’  electron  or  hole 
concentration.  This  localization,  however,  clearly  goes  beyond  impurity  scattering  of 
Dirac  fermions  considered  in  current  quantum  transport  theory10,12’25.  In  any  case,  a 
complete  theory  must  both  explain  the  reduction  of  Drude  weight  and  provide  a 
consistent  description  of  the  absorption  changes  in  interband  transitions. 

Our  infrared  spectroscopy  shows  that  the  frequency-dependent  response  of  Dirac 
fermion  in  graphene  can  be  described  by  the  Drude  form,  which  becomes  very  strong  at 
longer  wavelengths.  In  the  THz  range  ( /L  ~  300 jum ),  a  gated  graphene  monolayer  can 
absorb  over  15%  of  the  incident  radiation,  suggesting  that  graphene  can  potentially  be  a 
useful  new  THz  material.  More  importantly,  the  independent  determination  of  Drude 
weight  and  scattering  rate  for  the  free  carriers  in  graphene  has  allowed  us  to  discover  that 
the  free  carrier  Drude  weight  is  significantly  smaller  than  the  widely-accepted  theoretical 
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predictions,  but  the  sum  rule  for  integrated  intraband  and  interband  oscillator  strength  is 
strictly  obeyed. 

Methods 

Following  the  procedure  in  Ref.  20,  we  grew  graphene  on  copper  films  using  CH4  as 
the  feed  gas,  which  was  then  transferred  with  PMMA  support  to  a  Si/SiC>2  wafer  after 
wet-etching  to  remove  the  copper  film  by  FeCl3.  After  dissolving  the  PMMA  in  acetone 
solution,  high-quality  graphene  of  cm  size  on  the  Si/SiC>2  wafer  was  obtained.  (Fig. la) 
Raman  spectra  show  that  the  sample  was  indeed  monolayer  graphene.  Subsequently, 
Au/Ti  electrodes  (thickness  ~  50  nm)  were  deposited  in  vacuum  through  stencil  masks 
onto  the  graphene  sample  for  electrical  measurements.  The  doped  Si  substrate  (p-type, 
resistivity  ~  10  Q -cm)  underneath  the  290  nm  SiC>2  served  as  a  back  gate,  which  allowed 
us  to  tune  the  charge  carrier  density  in  the  graphene  sample.  All  optical  and  electrical 
measurements  were  performed  in  vacuum  (~  0.1  mTorr)  at  100  K. 
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Figures: 


Figure  1  Properties  of  CVD-grown  graphene  device,  a,  Optical  microscope  image  of 
CVD-grown  graphene  transferred  to  a  Si02/Si  substrate,  b,  Graphene  DC  conductivity  as 
a  function  of  gate  voltage.  The  conductivity  minimum  at  Vg  =  14V  defines  the  charge 
neutral  point  (CNP).  c,  Gate-induced  change  of  IR  transmittance  -AT/T  through  graphene 
at  Vg  =  -70  V  (hole  doped)  compared  to  transmittance  at  the  CNP.  The  spectrum  shows 
an  increase  of  free-carrier  absorption  at  low  wavenumbers  and  a  reduction  of  interband 
absorption  at  higher  wavenumbers,  d,  An  illustration  of  intraband  (i.e.  free-carrier 
absorption,  purple  arrow)  and  interband  (blue  arrow)  transitions  in  hole-doped  graphene. 
Intraband  absorption  increase  with  carrier  doping,  while  interband  transitions  up  to  2EF 
become  forbidden  due  to  empty  initial  states,  as  observed  in  c. 

Figure  2  Free-carrier  responses  in  graphene.  a,b,  Gate-induced  change  of  AC 
conductivity  Act'  in  hole-  and  electron-doped  graphene,  respectively,  for  30  cm_1<  to 
<450  cm'1  (solid  lines).  The  frequency  dependence  of  AC  conductivity  can  be  fit  by  the 
Drude  model  (dashed  lines),  c,  The  scattering  rate  T  and  ,d,  Drude  weight  D  at  different 
electron  and  hole  concentrations.  Surprisingly,  the  Drude  weight  is  substantially  lower 
than  theoretical  predictions  based  on  the  Boltzmann  transport  theory  (blue  line  in  d)  and 
shows  an  electron-hole  asymmetry,  e,  DC  conductivity  in  the  form  of  D/tiT  obtained 
from  the  Drude  fit  of  terahertz/far-IR  measurements  (red  dots)  agrees  well  with  the 
values  from  direct  DC  transport  measurements  (black  dots). 


13 


Figure  3  Interband  transitions  in  graphene.  a,b,  Gate-induced  change  of  interband 
AC  conductivity  Act'  in  hole-  and  electron-doped  graphene,  respectively,  over  600  cm 1 
<  oo  <  6000  cm1.  The  AC  conductivity  decreases  due  to  empty  initial  states  (hole  doping) 
or  filled  final  states  (electron  doping)  for  interband  transitions  below  2EF.  (The  sharp 
features  around  1200  cm'1  are  extrinsic  and  arise  from  SiC>2  phonon  resonances  of  the 
substrate).  Observed  interband  conductivity  decrease  is  appreciably  less  than  the 
theoretical  predicted  value  of  ne 2  /  2 h  (horizontal  dashed  line),  and  shows  electron-hole 
asymmetry,  c,  Upon  charge  carrier  doping,  integrated  conductivity  change  has  an 
increase  from  intraband  transitions  (red  symbols)  that  equals  the  decrease  from  interband 
transitions  (black  symbols).  The  two  values  agree  quantitatively  at  all  gate  voltages,  as 
required  by  the  oscillator  strength  sum  rule. 
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Supplementary  Information: 

The  frequency-dependent  conductivity  change  of  graphene  Act  at  different  gate 
voltages  is  obtained  from  the  measured  transmission  spectra  using  a  perturbation 
treatment,  because  a  monolayer  graphene  only  absorbs  a  small  fraction  of  the  light. 
Within  this  approximation,  the  complex  AC  conductivity  change  is  related  to  the 
difference  transmission  spectra  by 
-AT  4  n 

- (o))  =  —  Re(Acr  x  L)  (S_Eq-  1) 

T  c 

where  -AT IT  is  the  normalized  transmission  difference,  L  the  local  field  factor,  and 
Act  the  gate-induced  conductivity  change  in  graphene. 

For  suspended  graphene,  the  local  field  factor  is  one  and  we  obtain  the  well  know 
-AT  4  n 

form  of - (of}  = - Acr'(Sl).  In  our  device,  the  graphene  is  sitting  on  a  SiOi/Si 

T  c 

substrate  and  the  local  field  factor  L  is  uniquely  defined  by  the  refractive  index  of  Si  and 
SiC>2,  as  well  as  the  thickness  of  SiC>2,  and  it  can  be  determined  accurately  without  any 
adjusting  parameters.  Specifically,  L  has  the  form 


2  ikn„d 


L=l+r 


ag 


+ 1  T  t 

ag  gs  ga 


2iknad 

-  p  Syr 

'  gsga 


(S_Eq.  2) 


where  r  ,  r  and  r  are  the  Fresnel  reflection  coefficients  at  air-silica,  silica-silicon  and 

silica-air  interferences,  t  and  tga  Q^re  the  corresponding  transmission  coefficients,  ng 

and  d  are  complex  reflective  index  and  thickness  of  the  SiC>2  layer,  respectively,  and  k 
is  the  wavevector  of  the  incident  light. 

The  real  and  imaginary  part  of  AC  conductivity,  a'  and  a" ,  is  connected  by 
Kramers-Kronig  (K-K)  relation 
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\  ~zco  f 
cr  (®)= - 

7 T  J 


-2co  r  o'(co)dco 


TC  *  CO1  -CO2 


(S_Eq.  3) 


Because  the  K-K  relation  holds  true  for  all  the  gate  voltages,  the  gate-induced 
conductivity  change  Ao=Ao'+\Ao"  also  satisfies  the  same  K-K  relation.  Making  use  of 
the  S_Eq.  1  and  the  K-K  relation,  we  obtain  both  the  real  and  imaginary  part  of  the  gate- 
induced  AC  conductivity  from  infrared  transmission  spectra. 

The  gate-induced  infrared  absorption  from  the  p-doped  Si  substrate  is  negligible.  Real 

2  t-i 

n-e  r 


part  of  hole  AC  conductivity  in  silicon  is  described  by  o'  =  ■ 


mh  T2  +  co 2 


.  For  the  same 


carrier  density,  it  is  more  than  one  order  of  magnitude  smaller  than  graphene  conductivity 
across  the  experimental  spectral  range  (30  cm'1  <  co  <  6000  cm"1)  due  to  the  large  hole 


effective  mass  (  mh=  0.36  m0  )  and  small  scattering  rate  (T  ~  4  cm'1  at  100  K)  (S2). 
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